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Abstract 

The coupled system of renormalized Dyson- Schwinger equations for the 
electron self-energy and the photon propagator are supplied with the tree 
level vertex as Ansatz for the renormalized three point function. The system 
is investigated numerically. In the case of a massive electron, the theory 
is "weakly renormalizable" , i.e. cutoff independent for values of the cutoff 
below an upper limit. In this regime of cutoff independence, the quenched 
approximation yields good results for the electron self-energy. In the chiral 
limit, a logarithmic cutoff dependence of the electron self-energy is found. 
The question, whether a regime of cutoff independence with a spontaneously 
broken chiral symmetry exists in strongly coupled QED, remains open. 
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1 Introduction 



The occurrence of spontaneous chiral symmetry breaking in strongly coupled QED is 
one of the most challenging issues of non-perturbative quantum field theory nowa- 
days. In his pioneering work, Miranski reported a regime of cutoff independence 
(CI-regime)[| corresponding to a phase of QED with spontaneous broken chiral sym- 
metry Hfl. His results are based on a study of the Dyson- Schwinger equation for 
the electron self-energy, where electron loop corrections to the photon propagator 
are neglected (quenched approximation). In the quenched approximation, no renor- 
malization of the electric charge is required, and spontaneous symmetry breakdown 
occurs, if the bare charge exceeds a critical value |Q 0, ||. Subsequently, it was ar- 
gued by Bardeen, Leung and Love that a Nambu-Jona-Lasinio type interaction ||, 
must be included in order to render the renormalization of the Dyson-Schwinger 
equation consistent with the known renormalization properties of the theory ||. 
The results of the quenched approximation are caught into question, since vacuum 
polarization effects might enforce the renormalized coupling to vanish in the the 
infinite cutoff limit B (triviality). First investigations of this effect were done by 
Kondo by parametrizing the effect of the vacuum polarization §. Further insight 
in the issue of triviality was gained by the important work of Rakow . He numer- 
ically solved the coupled set of (bare) Dyson-Schwinger equations for the electron 
self-energy and the photon propagator, and therefore included vacuum polarization 
effects self-consistently. He finds a second order chiral phase transition and zero 
renormalized charge at the critical point J7|. This seems to rule out the quenched 
approximation, since the quenched approximation predicts an interacting infinite 
cutoff limit. 

Parallel to the studies of the truncated Dyson-Schwinger equations, extensive lattice 
investigations were performed in order to clarify the triviality problem of QED || 



|9], [T(| . Lattice simulations are not bounded by any approximation, but suffer from 
a small correlation length compared with the intrinsic fermionic energy scale due 
to the finite lattice sizes. In order to overcome the difficulty of a small correlation 
length, one fits the lattice data to two Ansatze of the equation of state (the bare 
electron mass as function of the fermionic condensate). One of these Ansatze favors 
an interacting infinite cutoff limit [[J , whereas the other Ansatz is compatible with 
triviality |K| . The present status is that the lattice data are not conclusive enough 



in order to distinguish between the two equations of state ||. 

In this paper, we further develop the results of and study the coupled set of 
renormalized Dyson-Schwinger equations for the electron self-energy and the pho- 
ton propagator. The tower of Dyson-Schwinger equations is truncated by using the 
tree level electron-photon-vertex as Ansatz for the renormalized three point func- 

1 We do not use the notion of a "scaling limit" in order to avoid confusion with the term in solid 
state physics, where it corresponds to the case of exact scaling (zero masses). 
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tion. In order to introduce our notation and to make the renormalization procedure 
transparent, we first derive the renormalized Dyson-Schwinger equations (section 0). 
All Ansatze are made for renormalized quantities. This will later turn out to be cru- 
cial in order to compare the full results with those of the quenched approximation. 
We first concentrate on the case of a massive electron (section |3|) . A cutoff indepen- 
dence is found for sufficiently small values of the cutoff. The theory is called to be 
weakly renormalizable. In contrast, the quenched approximation always predicts a 
Cl-regime with no upper limit on the cutoff. Although the quenched approximation 
is incapable to predict the correct Cl-behavior, it provides good results for physical 
quantities in the Cl-regime (subsection |3.2|) . We then focus onto the chiral limit 
(section ^). A logarithmic cutoff dependence is seen as well in the electron self- 
energy as in the vacuum polarization. The problem of the existence of a Cl-regime 
characterized by a spontaneous breakdown of chiral symmetry is still unsolved. 



2 Renormalized Dyson-Schwinger equations 

The generating functional for bare Green's functions of QED in Euclidean space is 
given by the functional integral 



W\j*,r] B ,V B }(eB,m B ) = J VA^VqVq expj- J d A x 



Ln(x) (1) 



- r] B (x)q(x) - q(x)r] B (x) - j B (x)A^(x) 

[A]{x)F ia ,[A]{x) + q(x)(i$ + irriB)q(x) + q(x)4-(x)q(x) , (2) 

where F Mi/ LA](a;) = d fl A v (x) — d v A^(x) is the field strength tensor, and e B , mg 
represent the bare electric charge and the bare electron mass respectively. At zero 
external sources, the generating functional W[0, 0, 0}(m B , e B ) is invariant under U(l) 
gauge transformations of the integration variables, i.e. 

q(x) — ► exp{ia(x)} q(x) , A^x) -> A^x) + d^x) . (3) 

Bare Green's functions are divergent implying the theory must be regularized and 
renormalized in order to make sense. Renormalization is performed by absorbing 
the divergences into the renormalization constants i?2,3,e,mj which relate the bare 
sources and the bare parameters to the renormalized ones, i.e. 

rj B (x) = Z- 1/2 rj(x), fj B (x) = Z~ l/2 fj(x) , j B (x) = Z~ 1/2 j,(x) , (4) 
e B = Z e e R , m B = Z m m R . 
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The generating functional for renormalized Green's functions is obtained from 
, r] B , f] B }(eB, m B ) by replacing the bare sources and bare parameters by the 
renormalized sources and the renormalized parameters respectively, i.e. 



W R [j^ rj, rj](e R , m R ) = W[Z 3 1/2 Z e 1 j fll Z 2 l/ \ Z 2 1/z r]}(Z e e R , Z m m R ) . (5) 

The renormalized Green's functions are obtained by taking functional derivatives 
of W R (H) with respect to the corresponding external sources. For example, the 
renormalized electron propagator S R and the renormalized photon propagator D B U 
are given by 



V2„ 



-V2=i 



St](x) 5r](y) 



j,v>v=o 



(6) 



After a change of the integration variables, the generating functional (|5|) can be cast 
into 



W R [j^r],T]](e R ,m R ) 



VA„ Vq Vq exp 



d 4 x 



L(x) 



- rj(x)q(x) - q(x)r](x) - j^xjA^x] 



(7) 



L(x) 



4e| 



F^A^F^A^x) + Z 2 q(x)iQq(x) + Z im R q(x)q( 



x 



+ Z 1 q(x)4(x)q(x) 
where we have introduced 



ry ry 

£>2 ^3 



Zi Z m . 



(9) 



(10) 



The generating functional T for renormalized one-particle irreducible functions (in 
the following vertex functions) is obtained from W R ([5]) by a Legendre transforma- 
tion, i.e. 



T[A^,ip,ip} = lnWfl[7 M , 77, 77] - j d A x j^A^x) + r](x)^(x) + ij;(x)r}(x) 
where the external sources are implicitly related to the fields A, ip, V> by 



5 In W R [j^ 7], rj] = 5\nW R [j fl: rj,r]] ^ = 5 In W R \j^ rj, rj\ 



Sr](x) 



5r](x) 



(12) 



Below we will use the renormalized electron photon vertex, which is defined by 



Exploiting the fact that the functional integral WrU^, tj, rf\ (|7|) is not changed by 
a shift of the integration variables A^, q, q generates the renormalized Dyson- 



Schwinger equations |TTJ. In particular, one finds 
(S R (x,y))~ l = [Z 2 i$ + Z im R ] 5(x - y) 



Z x l d A z d 4 "w^ pL D^ v (x,z)S R {x 1 w)A u {z,w 1 y) 



Z, 



6p,6(x -y) = —{-d5 m + d^d a ) D au (x, y) 



(14) 



(15) 



+ Z\ I d A z d A w d A v tr {^^S R {x, z)A a (v, z, w)S R (w } x)} D^ u (v, y) . 



Throughout this paper, we will truncate the tower of Dyson- Schwinger equations 
by making an Ansatz for the renormalized vertex function A u . In this case, the 
renormalized electron propagator and the renormalized photon propagator can be 



obtained by solving the coupled system (p^JT5| ). 

In order to investigate the consistency of the Ansatz, one studies the Ward-Takahashi 
identities [O, which provide relations among Green's functions induced by gauge 
invariance. One first realizes that Z\ = Z 2 must hold in the renormalized Lagrangian 
(|]) in order to preserve the invariance under the transformation (^). Combining 
(|1C|) with (f|), one finds that the electric charge is renormalized by the photon wave 

1/2 

function renormalization constant, i.e. e# = Z 3 e^. If the bare charge e# acquires a 
finite value in the limit of large cutoff (as suggested by the quenched approximation), 
and if further this limit enforces Z 3 to vanish, then the theory is only consistent with 
a zero renormalized electric charge. This scenario is referred to as triviality of QED 
in the literature. 

Exploring the fact that the generating functional Wr (0) is invariant under small 
gauge rotations (^), one obtains 



d { a ] A^(x,y,z) = iS R 1 (y,x)5(x - z) - iS R 1 (x,z)5(x - y) . 



(16) 



Once the coupled system ( TT4] JT5| ) was solved for a particular choice of the renor- 
malized vertex A^, one inserts the solution for Sr into (|16D in order to check the 
accuracy of the Ansatz for the vertex function A M . 



In the following, we will work in Landau gauge and set Z\ = 
invariance. We will study the Ansatz 

A^(z,x,y) = 7 M 5{z - x) S(z - y) 



Z 2 as imposed by gauge 



(17) 



for the renormalized vertex function, which is the tree level electron-photon vertex. 
Recently, the full one-loop QED vertex was obtained |]12| . The results might provide 
the structure for a more general ansatz than (|T7|). Note that the choice of the 
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renormalized vertex is the only freedom we have. There is no further possibility to 
argue in favor of a four fermion interaction first introduced in 0. Note also that 
a four fermion interaction arises and can be naturally incorporated in renormalized 
Dyson-Schwinger equations in the dual formulation of QED ||13[| . 
The coupled Dyson-Schwinger equations are then solved by parametrizing the renor- 
malized electron propagator and the renormalized photon propagator in momentum 
space by 

c n 1 n-R f \ 4nD R {p 2 ) ( p^pA 

Up to the occurrence of the renormalization constants Zi^, the derivation yields 
the same results as reported in |7|]. The first Dyson-Schwinger equation is straight- 
forwardly reduced to two equations to determine F(p 2 ) and S(p 2 ). Using a sharp 
0(4)-invariant cutoff A to regularize the momentum integration, these equations are 

S(p 2 ) = m + Z X f dk 2 k 2 [" d9 sin 2 9 3§ , (19) 



F(p 2 ) = Za 



1 



TT 2 p 



(20) 



x 



A ^fc 4 f d0 sin 2 9 F ^J ° R [ q2) hq 2 cos 9 - 2fcp sin 2 9} 
o jo s(/r) g 4 L > 



where we have introduced tuq = Z m mR, s(k 2 ) = F(k 2 )k 2 + Y,(k 2 ) and q 2 = k 2 + p 2 — 
2kpcos9. The equation for the photon propagator fll5]) needs more thoughts. The 
vacuum polarization 11^ (p) due to the electron loop is transverse, i.e. p^ li^ip) = 0, 
if the regularization prescription does not violate gauge invariance. Unfortunately, 
the sharp 0(4)-invariant cutoff prescription spoils gauge invariance implying that 
the vacuum polarization IT^Qo) acquires spurious terms proportional to 5^, which 
are, in addition, quadratic divergent. On the other hand, regularization schemes 
consistent with gauge invariance are very time consuming when solving the cou- 
pled equations (|1^J15|) numerically. In order to circumvent this problem, one starts 
calculating Xl^ip) using a regularization scheme which respects gauge invariance 
(e.g. Schwinger proper time, Pauli-Villars regularization). The result for U^, (jp) is 
the transverse projector times a function depending on p 2 . In this function, one first 
introduces a second regularization which corresponds to the 0(4)-invariant cutoff 
scheme, and then removes the regulator of the scheme which is compatible with 
gauge invariance. The final result allows a fast numerical treatment and is consis- 
tent with the transversality of the vacuum polarization {jp). In the limit of a 
large cutoff, all regularization schemes which respect the symmetries of the model 
are supposed to yield the same result for physical quantities. In agreement with the 
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result stated in 0, this lengthy procedure yields 

Dr{p 2 ) a R 3n 2 

x f dk 2 k 2 r d9 sin 2 6 F [ q t] F{q ^ |^(8cos 2 g-2)-3/2| , 
Jo Jo S {Q+) S {Q-) [P J 

where a R = Cr/Att and q± = (k ±p/2) 2 . 

3 QED with massive electrons 

In order to calculate the renormalization constants and mo as function of the 
cutoff, we impose renormalization conditions. In order to fix Z\, we demand that 
the inverse renormalized electron propagator has the canonical kinetic term, i.e. 
F(0) = 1. Our conventions imply that the residue of the renormalized photon 
propagator D R (0) can be identified with a R . Providing a value for a R determines 
the constant Z 3 . Finally, we set the scale of the electron self energy S(p 2 ) by the 
constraint S(0) = So, which yields a condition to calculate ttiq. Although mo might 
tend to zero for infinite cutoff, the renormalized (current) mass m R is non-zero. 
The behavior of m Q in the quenched approximation might serve as an illustrative 
example ||, ||], i.e. 



/ u\ w / 3e 2 

m ° = \aJ mR ^ ' with W = 1 ~ ] J 1 ~4^' ( 22 ) 

where /i is an arbitrary renormalization point. Since mo will turn out to be different 
from zero for any fixed value of the cutoff, the renormalized (current) mass m R is 
also non-zero implying that the above set of renormalization conditions corresponds 
to the case of a massive electron. The chiral symmetry is explicitly broken. 

3.1 Numerical results 

The set of coupled equations (p!9| , pO| , pTD for the functions F(p 2 ), S(p 2 ), D R (p 2 ) was 



numerically solved by iteration. The integrals were done by Simpson's integration 
with a step-size control. A Cl-regime is found for sufficiently small cutoff. A repre- 
sentative result is shown in figure 1 for a R = D R (0) = 0.35. The function F(p 2 ) is 
weakly p-dependent. The maximum deviation from one occurs for p 2 = A 2 and is 
0.04776 for InA 2 /!! 2 , = 16. The self-energy roughly stays constant for hip 2 /!! 2 , ~ 4 
and smoothly decays afterwards. The vacuum polarization also stays constant up 
to approximately the same momentum and then decays according to 

— = CiK) - c 2 (a R ) \nA 2 /p 2 , for (lnp 2 /£ 2 > 4) , (23) 
D R {p 2 ) 
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where ci(0.35) ~ 2.97 and where Ci ~ 1/tt 2 for a wide range of renormalized cou- 
plings. The renormalization constants are given in the table below. 



lnA 2 /£g 


Zi 


^3 


m /S 


7 


1.0097 


0.75537 


0.49437 


10 


1.0183 


0.64295 


0.33425 


13 


1.0305 


0.52927 


0.20643 


16 


1.0497 


0.41307 


0.10936 



The renormalization group flow is shown in figure |2] for o.r = 0.35. The vertex 
renormalization constant Z\ (which is identical to the wave function renormaliza- 
tion constant Z 2 of the electron) stays close to one. The photon wave function 
renormalization constant almost decays logarithmically. The behavior of the bare 
mass tuq corresponds to a powerlaw decay with logarithmic corrections. 
In order to check the accuracy of the Ansatz ([17]) for the vertex, we study the Ward 
identity (|i~6"D in momentum space and at zero photon momentum, i.e. 



A„(0,p) = |-S 

A* 



-1 



(P)- 



(24) 



Using the parametrization (|1|) of the electron propagator and neglecting derivatives 
of the functions F(p 2 ) and S(p 2 ), the Ward identity (B4T) becomes 



A M (0,p) = F{ P 2 ) lll 



(25) 



The deviation of F(p 2 ) from unity measures the violation of the ward identity due 
to the Ansatz (O) for the vertex. The numerical result shows that this violation is 
small, since the deviation of F(p 2 ) from unity is always beyond 5% for the values of 
the cutoff used in table 1. 

Increasing the cutoff for fixed renormalized coupling cur, one observes a critical upper 
limit Ac of the cutoff. Beyond this critical value, no solution of the set of equations 
(|l^,[l^,^(]) was found. This result can be anticipated from the asymptotic behavior 
of the vacuum polarization (|2"3"D. For a fixed (positive) value of the renormalized 
coupling cur = Dr(0), equation (p3| ) cannot be satisfied for an arbitrarily large cutoff 
A. Figure 3 shows the relation between the maximal possible renormalized coupling 
and the critical cutoff. The results indicate that in the limit A — > oo a solution of 
the coupled Dyson- Schwinger equations only exists for —* 0. This phenomenon 
is referred to as triviality of massive QED. Note, however, that for values of the 
cutoff smaller than the critical value the electron self-energy is cutoff independent. 
The theory is called weakly renormalizable. For moderate values of the renormalized 
coupling, the theory allows for a large value of the critical cutoff and therefore for 
a cutoff many times bigger than the typical energy scale set by the mass of the 
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electron. This implies that weakly renormalizable QED is perfectly compatible with 
the observations in nature. The lack of a cutoff independence at high momentum 
does no harm, since QED in the real world is embedded in the Weinberg-Salam 
model at high energies. 



3.2 The quenched approximation 



Neglecting fermion loop effects in the photon propagator is called quenched approx- 
imation. In our formulation, this corresponds to replacing (|2l"D by 



D R (p* 



(26) 



The vacuum polarization Dr(p 2 ) is momentum independent, and from the renor- 
malization conditions (see beginning of section ^) one immediately obtains Z% = 1. 
For a constant vacuum polarization Dr, the angle integral in ( ^0|) can be done ex- 
plicitly and yields zero. This implies that F(p 2 ) = Z\ is also constant, and taking 
into account the renormalization conditions, we have Z\ = 1. The only remaining 
non-trivial equation is the integral equation for the electron self-energy flT9p. For 
standard QED (massless photon), this integral equation can be transformed into a 
non-linear differential equation with appropriate boundary conditions [[l], |], || f|. 
In the table below, we compare the renormalization constants of the full and the 
quenched approach for an = 0.35 and a cutoff in the Cl-region lnA 2 /£Q = 16. 







Zz 


m /% 


quenched 
full 


1 

1.0497 


1 

0.41307 


0.23172 
0.10936 



Large deviations are found. Note, however, that renormalization constants are not 
physical observables, and the quenched approximation might improve when physi- 
cal quantities are studied. Figure |] shows the electron self-energy in the quenched 
approximation in comparison with the full result. One finds that the quenched 
approximation yields good results for the electron self-energy at least for small mo- 
mentum p. This behavior can be understood by a scaling argument. Due to cutoff 
independence, the momentum dependence of the electron-self-energy for a small mo- 
mentum p and for cutoff slightly below the upper critical value is also obtained in 
a calculation with a cutoff far beyond the critical cutoff and momentum p compa- 
rable to this cutoff (see figure 1). At small cutoff, however, the bare electron mass 
m is large (see table 1), and the electron loop contributing to the inverse vacuum 
polarization in fl21~|) is negligible, hence the quenched approximation is good. 
If the cutoff exceeds the upper critical limit, a solution of the full equations ceases 
to exist, whereas the quenched approximation still predicts a solution. This implies 
that polarization effects are important to address the Cl-behavior of the model. 
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We conclude that within the quenched approximation one cannot decide whether 
a Cl-regime exists or not. If a Cl-regime exists, the quenched approximation may 
be a good approximation for physical observables, although non-observables (e.g. 
renormalization constants) might turn out to be completely different. 



4 Cl-violation in chiral symmetric QED 

In the following, we will study the chiral limit = (implying m = 0) of QED, 
which requires a new set of renormalization conditions. Again we demand that the 
residue of the electron propagator is one, i.e. F(0) = 1, which fixes Z\. We are 
interested in a phase with spontaneously broken chiral symmetry. We therefore 
insist on E(0) = So, which must now be accomplished by a choice of Z3 since mo 
is now identical zero. All renormalization constants are fixed. The renormalized 
coupling olr = -D_r(0) will be self-consistently calculated. 

We numerically solved the coupled equations ( |l9| , p0l , pl~|) for this set of renormal- 
ization conditions. The result for the electron self-energy X(p 2 ) and the vacuum 
polarization D^{p 2 ) is shown in figure 5, which should be compared with figure 1. 
We find a logarithmic dependence on the cutoff. Although the qualitative behavior 
of 1/Dji(p 2 ) is qualitatively the same as in the case of massive electron, the plateau 
at small p 2 does not stabilize, but continuously increases with increasing cutoff. 
One might think that this Cl-violation in the self-energy is due to the change of 
the plateau of Dr{j> 2 ) when the cutoff is varied. One might therefore be tempted to 
parametrize the vacuum polarization by 

' Cz{a R ) - c A {a R ) lnA 2 /p 2 , for A 2 -> 00 (27) 



D R (p 2 ) 

and to search for a Cl-regime by solving the reduced set of equations ( |I9| , |2"0"| ) . How- 
ever, it turns out that the Cl-violation in the electron self-energy qualitatively re- 
mains the same. 

Two possible explanations of this cutoff dependence are immediate: the theory is 
trivial in a sense that the limit A — > 00 is only compatible with olr = 0. In this case, 
the interaction is not strong enough to spontaneously break chiral symmetry, and 
one cannot keep the renormalization condition E(0) = S 7^ 0. A second possible 
explanation is that in the chiral limit the solutions are sensitive to the infrared 



behavior of the integrals in ([T{],g(j,£l]) . The momentum independent Ansatz for 
the renormalized vertex fllTD might be to crude for the chiral limit and induces the 
logarithmic Cl-violation. An improved Ansatz for the vertex function might provide 
a Cl-regime with a spontaneous broken symmetry. None of these two explanations 
can be favored without further studies beyond the Ansatz (|17|) for the vertex studied 
here. 
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5 Conclusions 



We numerically studied the coupled set of renormalized Dyson-Schwinger equations 
for the electron and the photon propagator using a tree level electron-photon- vertex 
as an Ansatz for the renormalized vertex function. We imposed the constraint 
Z 2 = Z\ between the electron wave function renormalization constant Z 2 and the 
vertex renormalization constant Z\ as required by gauge invariance. The violation 
of the Ward identity is smaller than 5%. 

In the case of a massive electron (explicitly broken chiral symmetry) and fixed 
renormalized coupling cur, a regime of cutoff independence (Cl-regime) was found, 
if the cutoff is below an upper critical value. The theory is weakly renormalizable. 
For larger values of the cutoff than the critical value, no solution of the Dyson- 
Schwinger equations exists. The results indicate that the infinite cutoff limit is 
only compatible with a zero renormalized charge. In the Cl-regime, the quenched 
approximation yields reasonable results for the electron self-energy, although the 
values of the renormalization constants completely differ from the full result. 
In the case of chiral symmetry (zero renormalized mass of the electron) a logarithmic 
dependence on the cutoff was found. The question, whether this cutoff dependence 
is due to the crude Ansatz for the renormalized vertex or induced by triviality, 
cannot be answered at the present stage of investigations. The question, whether 
a Cl-regime with spontaneously broken chiral symmetry exists in strongly coupled 
QED, cannot be answered from the viewpoint of the truncated Dyson-Schwinger 
equations supplied with tree level vertex. 

Acknowledgements: I thank J.-B. Zuber for helpful discussions. 
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Figure captions 



Figure 1: Cl-behavior: the electron self-energy £(p 2 ) and the vacuum polarization 
Dr as function of the momentum p 2 for several values of the cutoff A in units of 
S = S(0). 



Figure 2: The renormalization group flow for a R = 0.35; -Zi(A) short dashed line, 
Z 3 (A) long dashed line, m (A)/E solid line. 

Figure 3: The maximal possible renormalized coupling a R as function of the critical 
cutoff A c (left). The ratio Z 3 /an at the critical cutoff (right). 



Figure 4: The electron self-energy £(p 2 ) in the quenched approximation compared 
with the full result for a R = 0.35 and A 2 = 8.9 10 6 E . 

Figure 5: CI- violations: the electron self-energy S(p 2 ) and the vacuum polarization 
D R as function of the momentum p 2 for several values of the cutoff A in units of 
S = S(0) in the chiral limit m R = 0. 
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